Numerical study of t2 g orbital system with ferromagnetic polarization 
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Finite temperature orbital state in a ferromagnetic Mott insulator with triply-degenerate tig orbital is inves- 
tigated numerically. We employ the quantum Monte Carlo simulation with the loop algorithm. Indications 
for conventional staggered-type orbital order are not remarkable down to the lowest temperature to which the 
present simulation can get access. Physical parameters monitoring the off-diagonal orbital order, which is 
characterized by a linear combination of the [dy^d^dxy) orbital- wave functions with equal weights, are not 
conspicuous. It is found that a orbital gap-like behavior appears in the uniform orbital susceptibility. This 
is supported by a threshold behavior in the staggered correlation function in a calculation with the additional 
Ising-type interaction. Some rigorous remarks for the long-range orbital order are also presented. 

PACS numbers: 75.40.Mg. 75.10.-b, 75.30.Et 
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Orbital degree of freedom and its strong coupling with spin, 
charge and lattice bring about a variety of fascinating physi- 
cal phenomena in transition-metal oxidesji*^ Colossal mag- 
netoresistance effect observed in vicinity of charge/orbital or- 
dered state in perovskite manganites is a typical example. 
Classical aspects of the orbital, e.g. long range orbital orders 
accompanied with the Jahn-Teller distortion and their roles on 
the magnetic exchange interaction, have been fully examined 
and been almost settled nowadays. On the other hand, quan- 
tum aspects of the orbital are still unrevealed and are expected 
to open a window for a new research field in correlated elec- 
tron system. 

Perovskite titanate /?Ti03 (R: rare-earth metal ion) is one of 
the candidates as a material where quantum orbital physics is 
realized. A formal valence of Ti ion is 3+, and one 3d electron 
occupies one of the triply degenerate t2 g orbitals. In this paper, 
we denote, for simplicity, the three t2 g orbitals, {d yz , d a , dxy) 
by {d a , dp, dy), respectively. Because of a weak Jahn-Teller 
coupling and the novel symmetry of the t2 g orbital-wave func- 
tions, quantum aspects of the orbital are expected to become 
evident4i££ 

In particular, a ferromagnetic Mott insulator with the ?2 g or- 
bital degeneracy attracts attention as a simple orbital system 
where spin and charge degrees of freedom are quenchediiSjE 
One representative material is YT1O3, where the isotropic 
spin- wave dispersion relation is observed by the inelastic neu- 
tron scattering^ This result seems to contradict a long-range 
orbital order confirmed by several experiments ! n ' 12,13 ' 14 The 
t2 g orbital Hamiltonian in a ferromagnetic Mott insulator at 
zero temperature was examined in Ref . |H Through the analyt- 
ical calculations, the two kinds of long-range orbital ordered 
states^ were proposed. It was shown that these orders set in 
at only at zero temperature due to strong quantum fluctuation. 
However, the puzzling spin-orbital properties in YT1O3 still 
remains an open issue. Even in a pure theoretical view point, 
finite-temperature (T) orbital states with large quantum fluc- 
tuation are unclear. 

In this paper, we study the t^ g orbital system with fer- 
romagnetic polarization. Our purpose is to examine finite- 
temperature orbital state in the idealized t2 g orbital model by 
using an unbiased method. We employ the quantum Monte- 
Carlo (QMC) simulation with the loop algorithm. Down to 



at least around T ~ 0.3/ with the coupling constant J, in- 
dications of the conventional staggered-type order, the off- 
diagonal orbital order and the orbital dimerized state are not 
conspicuous. Instead, the orbital gap-like behaviors are found 
in the uniform orbital susceptibility below around 0.87. 

We start from the Hamiltonian for the fag) system with 
ferromagnetic polarization in an ideal Perovskite lattice. The 
generalized Hubbard Hamiltonian in a simple-cubic lattice 
with the triply-degenerate t2 g orbitals is given as 
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We define the electron annihilation operator d^ a with or- 
bital 4 = (a,j3,y), spin s = (T,4) at site i, and the elec- 
tron transfer integral t?j . A pair of the nearest neighboring 
(NN) sites is denoted by (ij) . The intra-orbital Coulomb in- 
teraction U, the inter-orbital one U', the Hund coupling K, 
and the pair-hopping I are introduced. In an atomic limit, 
there is the relations U = U' + 21 and K = I. We also in- 
troduce a number operator, = Y*s n i$s = Hs d J^ s d i^s' which 
has a conservation = 1 ■ When the electron transfer- 

integrals through the O 2p K orbitals are only taken into ac- 

count, the transfer integral is written as a simple form fa = 
t8^i(8 ai ^ + 8y^). A subscript l(=x,y,z) indicates a direc- 
tion connecting i and j. We define two "active" orbitals, a/ 
and b[, with a finite transfer-integral, and an "inactive" one, 
ci, with no transfer-integral. That is, (a x ,b x ,c x ) = (j3,y, a), 
(a y ,by,Cy) = (y,(X,P), and (a z ,b z ,c z ) = (<X,fi,y). It is conve- 
nient to introduce the pseudo-spin (PS) operator for the active 



orbitals T • = i £ 



with the Pauli matri- 



ces <7. The effective Hamiltonian for the fa g ) system which 
has strong on-site Coulomb interactions is obtained by the 
perturbational processes ! 7 ' 9 ' 16 The virtual intermediate states 
are classified by the irreducible representations for the fag) 2 
states, i.e. l A\, 1 72, E and 3 7i. The lowest energy state is 
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3 7i with the energy of U' — K. We assume spins are saturated 
in ferromagnetic phase where 3 7\ is only the relevant inter- 
mediate state. The Hamiltonian studied in the present paper 
is 
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where the exchange constant is / = t 2 /(U' — K). This 
Hamiltonian is rewritten by using the PS operator as Jff' = 
47 {T' • T'j + (tijcjijc, — 1 ) /4}, and is the same Hamilto- 
nian studied in Refs.[7|and[9l In Ref. the authors proposed 
the two orbital orders where the wave-functions are the linear 
combinations of the (d a ,dp,dy) orbitals with equal weight^ 
which are called the off-diagonal type orbital orders in this 
paper. 

Before showing the numerical results, we remark some rig- 
orous results for the orbital order. There is a novel symmetry 
in this Hamiltonian: in each plane perpendicular to a direc- 
tion /, an electron number with the inactive orbital c/, Yli n ic v 
is conservedi^ 7 ^ A symbol Yl\ implies a sum of sites in a 
plane being perpendicular to I. The Hamiltonian is invari- 
ant under the two-dimensional U(l) gauge transformation de- 
noted by U{9) = exp(— /0^n IC/ ). By following the general- 
ized Elitzur's theorem ^ 17 ' 18 the physical quantities which are 
not invariant under this transformation has a vanishing mean 
value. This is proven from a combination of a theorem about a 
gauge transformation in a reduced dimension and the Mermin- 
Wagner theorem. 17 ' 1 - By using the theorem, we prove that the 
long-range order for the following orbital operator does not 
appear at finite temperature, 
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for / = (x,y,z) with numerical constants C a and Q,. That is to 
say, (Oi) = at finite temperature. 

To analyze the Hamiltonian in finite temperature by using 
an unbiased method, we employ the QMC method with the 
loop algorithm! 20 ' 21 There is no negative-sign problem. Sim- 
ulations are performed on L 3 cubic lattices with the periodic- 
boundary condition. The Suzuki-Trotter decomposition with 
the Trotter number n is adopted. To check efficiency of the 
present QMC simulation, we calculate the auto-correlation 
time^ 3 - for the parameter L^ 3 T x . The auto-correlation time 
grows up monotonically with decreasing temperature, and 
reaches around 750MC steps at T/J = 0.3. Thus, we adopt, in 
the following simulation, 2x 10 3 MC steps for thermalization, 
and 10 4 MC steps for measurements. We estimate averages 
and errors for the physical quantities in 10 independent simu- 
lations. System size is taken to be L = 6 — 12 in a cubic lattice. 
Numerical data obtained in the Trotter numbers «=10, 20 and 
30 are extrapolated. 

First, we show results for the staggered-type orbital order. 
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FIG. 1: (a) Staggered correlation function M x for the orbital PS 
operator Tf , (b) Staggered orbital susceptibility % x . For comparison, 



^/X X T/L 3 is a l so plotted in (a). 



We introduce the staggered correlation function along direc- 
tion I, 
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and the staggered orbital susceptibility 
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where — 6nL?, x is an imaginary time, T^(i) is the z compo- 
nent of the PS operator at T, and (• • • ) implies the QMC aver- 
age. In the classical limit, M l should be identical to y}T jV> . 
Numerical results of M x are presented in Fig. [TJ a )- We have 
checked that three M 1 for / = {x,y,z) coincide with each other 
within numerical errors. In L — 6 case, M x shows a broad 
peak around T /J = 1.5. This peak is, however, smeared out 
with increasing L. All obtained values are less than 10% of 
the maximum value, i.e. one. For comparison, \J % X T jl? 
is also plotted in the same figure. In high temperatures, M x 
and 'x x T '/L 3 merge together for all L as expected. A devia- 
tion between the two is remarkable below T/J ~ 1 where the 
quantum effect starts to be effective. In both the two quanti- 
ties, there are no divergent behaviors at T = 0. As shown in 
Fig. 01b), numerical data of % x for all system size L are scaled 
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Qi = T.l=(x,y,z) T L defined by 
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FIG. 2: (a) Correlation function G@ for the off-diagonal orbital op- 
erator Qj. (b) Correlation function G^" 1 ' for the orbital PS operator 
7g, System size in (b) is taken to be L = 10. Results in G x ^> and 
G*(4) jj-g multiplied by five. 



by a single curve in a whole temperature range. This implies 
that the correlation function £y TT /(-l) i+ ^'(r/(T)r^(T / )) in 
Eq. © is of the order of L? . That is to say, this correlation 
is rather short ranged down to at least 0.37. The obtained % x 
remains to be finite down to the lowest temperature. 

This decrease in M 1 and % below T/J ^ 1 may be ex- 
plained by a scenario that the orbital order of the off-diagonal 
operator proposed in Ref. is developed in low temperatures; 
in both the two wave-functions claimed in Ref. |9|, the cor- 
relation function (T/.TjS) for every NN pairs of sites i and j 
vanishes. To examine this possibility, we calculate the physi- 
cal quantities for the off-diagonal orbital operator. In general, 
thermal average of the two-point function for an off-diagonal 
operator may be calculated by utilizing the improved estima- 
tor in the loop algorithm! 20 i 21 In the present QMC simula- 
tion, however, this method is not realistic. This is because 
a value of the improved estimator depends on a shape of the 
loop, since the orbital interaction explicitly depends on the 
bond direction. Instead, the NN correlation functions for the 
off-diagonal operators are able to be calculated, when this 
is included in the Hamiltonian explicitly. We calculate the 
following correlation function for the off-diagonal operator 
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where z = 6. For the orbital ordered states proposed in Ref. |9l 
G^ = 0.22 and 0.16 for the type-I and type-II orders. Nu- 
merical results of G^ presented in Fig. |2{a) show weak size 
dependence and do not show remarkable increasing with de- 
creasing temperature. The obtained values are about 35% of 
the values expected from the ideal off-diagonal orbital orders. 

Further information for the off-diagonal orbital order is ob- 
tained by the staggered correlation as a function of distance: 
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where £/y\ represents a sum for the m-th NN sites i and j 

along direction /. We calculate G'^'s for / = x up to m = 4 
[see Fig. EJb)] . With decreasing temperature, G*^ 1 ' and G x ^ 
monotonically increase even below T/J ~ 1 where M and % l 
in Fig. Q] start to decrease. On the other hand, reductions are 
seen in G x ^ and G^ 4) around T/J = 0.8. That is, the de- 
creases in M' and X are attributed to those in the long-range 
correlations of TL and short-range correlations still remain 
to grow up. The obtained result also suggests that the de- 
creases in M l and X are n °t attributed to development of the 
off-diagonal orbital correlation; when the correlation for the 
off-diagonal operator Qj grows up at low temperature, reduc- 
tion in G^" ! ' should be remarkable in short range. That is, 
G'O and G l{ ~ 2 \ instead of G 1 ^ and G l( - 4 \ are expected to be 
reduced. This expectation is in contrast to the numerical re- 
sults in Fig. Efb). 

Development of the short-range orbital correlation shown 
above suggests a kind of the valence-bond state in low- 
dimensional quantum magnets. 24 However, we are able to ex- 
clude a simple orbital dimer state. We calculate the dimer 
order parameter defined as 
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with a unit vector e/ along direction /. Results of D 1 show 
monotonic decreases with decreasing temperature and in- 
creasing the system size (not presented in figure). 

The following results indicate that the low-temperature or- 
bital state shows a gap-like feature. The calculated uniform 
orbital susceptibility defined by 



Xu 



4 

NT \ 



({EE^)} : 



(9) 



[see in Fig.^a)] starts to decrease around T/J = 0.7 and tends 
to vanish at low temperatures. Size dependence is seen below 
T/J ~ 0.3, but a global gap-like feature in %[, is not sensi- 
tive to the system size. The temperature, where %[, starts to 
decrease, i.e. T ~ 0.77, almost coincides to the temperature 
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FIG. 3: (a) Uniform susceptibility for the orbital PS operator Tf . 
(b) Correlation function M z calculated in the Hamiltonian jtf' + Jtf' . 
System size and the Trotter number in (b) are chosen to be L = 8 and 
n = 10, respectively. 



where J' is a positive coupling constant. This term promotes 
the staggered-type orbital order characterized by M z . Results 
are shown in Fig.|5Jb) for several values of J' /J. A threshold 
value for J' /J seems to exist; an increase of M~ in low tem- 
peratures is not seen for J' /J = 0.05 and 0.15. This is con- 
sistent with finite values of % x at low temperatures shown in 
Fig. QIb). These results indicate a possibility that a quantum 
gapped state is broken by the Ising-type interaction as well 
known in the quantum critical issue. 



In summary, we present a numerical study of the finite- 
temperature t2 g orbital state in a ferromagnetic Mott insula- 
tor. Remarkable developments are not seen in the staggered 
orbital correlation and the NN correlation for the off-diagonal 
orbital operator at least down to T/J ~ 0.3. Short-range cor- 
relation remains to grow up in low temperatures, but a pos- 
sibility of the dimerized state is excluded. It is found that 
the uniform orbital susceptibility shows a gap-like feature be- 
low T I J ~ 0.7. This result is supported from a calculation 
in the model where the Ising-like interaction is added. One 
possible scenario for the low temperature orbital state is a 
short-range singlet state where singlet pairs are randomly dis- 
tributed in a crystal, and/or rearrangement of pairs occurs dy- 
namically. A singlet wave function between sites ; and j are 
(\diafdji,,) — \dib,dj a ,)) / v2 where concerning orbitals depend 
on the bond direction /, i.e. a directional singlet pair. To 
clarify a low temperature orbital state in more detail, QMC 
simulations in large cluster sizes with efficient algorithms are 
necessary. A variational-type approach based on the above 
directional singlet assumption may be also helpful. 



where the staggered susceptibility % x and the short-range cor- 
relation function G x ^ and G x ^ have broad peaks, as shown 
in Figs. 02b) and|2jb). This gap-like feature is also confirmed 
by the following calculation. We add the Ising-like interaction 
term into the Hamiltonian Eq. ©I 

Jf = 2/ £ Ut^ - n ia n ja - n ip njp) , (10) 
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